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Abstract 


and 


In the present paper, the shrinkage estimators have been 
developed based on the Upper record values under the 
Bayesian Minimax Risk Criteria. The properties of the 
proposed shrinkage estimators have been studied in terms of 
the relative efficiency with respect to improved estimator. 
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Introduction 

The Rayleigh distribution is a suitable model for life 
testing experiments and clinical studies. A significant 
amount of work has been done related to Rayleigh 
model. The origin and other aspects of this 
distribution are found in Siddiqui (1962) and Hirano 
(1986). Polovko (1968), Dyer & Whisenand (1973) 
demonstrated the importance of this distribution in 
electro vacuum devices and communication 
engineering respectively. Its many applications have 
also been discussed by Ariyawansa and Templeton 
(1984). Hendi et al. (2007) presented a Bayesian 
analysis for record values from the Rayleigh model. 
Recently, Prakash and Prasad (2010) have present a 
Bayes Prediction interval for the Rayleigh model. 

The objective of the present paper is to study the 
properties of the Shrinkage estimator based on the 
Upper record values under the Bayesian Minimax Risk 
Criteria. 

The Considered Model and Prior 
Distributions 

The probability density function and distribution 
function of the Rayleigh model with the parameter 0 
are given as 


F (x ; 0) = 1 - exp 


' x 2 ^ 


( 2 . 2 ) 


Here 0 is known as the shape parameter. The 
Reliability function 'F ( t ) and Failure rate p ( t ) 
for any specific mission time t > 0 are obtained 
respectively as 


and 


'P ( t ) = exp 


7 

v“9y 


; 0 > 0 


( 2 . 3 ) 


p{ t ) = ~;0>O- (2.4) 

0 

Here, we are not going in to the debate or to justify the 
questions of the proper choice of the prior distribution. 
For situations where life tester has no prior 
information about the parameter, they may use the 
quasi-density prior. Here we consider a quasi prior of 
the following form: 

71 ( 0 ) gc — !-j- ;d > 0, 0 > 0 . ( 2 . 5 ) 

0 


If d = 1 we get a non-informative prior (Jeffrey's 
(1961)) 

71 (0 ) oc — 

' 0 


and if d = 3 we get the asymptotically invariant prior, 
proposed by Hartigan (1964) 

7^ 2 ( 0 ) °c • 


Record Values and the Estimation 


f (x ; 0) 


2 x 

— exp 
0 
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; x > 0, 0 > 0 


( 2 . 1 ) 


Record values found in many conditions of daily life 
as well as in numerous statistical applications. The 
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record values, can be viewed as order statistic from a 
sample whose size is determined by the values and the 
order of occurrence of the observations. In a little over 
thirty years, a large number of publications devoted to 
records have appeared. 

This is possibly due to the fact that we encounter this 
notion frequently in daily life, especially in singling 
out record values from a set of other ones and in 
registering and recalling record values. 

Motivated by extreme weather conditions. Chandler 
(1952) introduced the study of record values and 
documented many of the basic properties of records. 
Feller (1966) gave some examples of record values 
with respect to gambling problems. Resnick (1973) 
discussed the asymptotic theory of records. The theory 
of record values and its distributional properties have 
been extensively studied by Ahsanullah (1995), Arnold 
et al. (1998). 

Let X j, X 2 ,„. be the sequence of independent and 
identically distributed random variables following the 
model (2.1). The observation X ■ is an upper record 
value if it is larger than all preceding observations, i.e., 
if X j > X j , for all i < j . The indices at which these 

upper record values occur are called upper record 
times { U (m), m > 0 }, where U (0) = 1 and U (m) 
= min {j : j > U(m - 1), X , >X U(m _ 1) }. 

Let X = X U(|) ,X U(2) ,...,X U(n) be the n upper record 

values from model (2.1). The likelihood function is 
thus obtained as 

L(xie)=f (x U(n) ,e)n />(x U(i) ) (3.i) 

i=l 

It follows from (2.1), (2.4) and (3.1) that 
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e>o,s = x 2 U(n) . (3.2) 


The Maximum likelihood estimator (MLE) and 
unbiased estimator for the parameter 0 based on the 
upper records values are given respectively as 

0ML =-and0 u =S. 

n 


Bayes Estimation 


Under the Bayes theorem, the posterior distribution of 
the parameter 0 is obtained as 


g (eix)= L MHe) 
8 ' ; jL(xie).u(e)de 


g n + d-l 

r (n + d - 1) 0 n+d 



-1 

ej 


;0 > 0, d > 0. 


The Bayes estimator for the parameter 0 under the 
SELF (squared error loss function) is simply the 
posterior mean and given as 


0 bs = © i S ; CO 


1 

n + d - 2 


(4.1) 


Based on the Bayesian context, it is standard practice 
for most statisticians to consider the squared error loss 
function as it does not lead to extensive numerical 
computation. The symmetric nature of this function 
confers equal importance to overestimation as well as 
underestimation. If the SELF is taken as a measure of 
inaccuracy then the resulting risk is often too sensitive 
to the assumptions about the behaviour of the tail of 
the probability distribution. In addition, in some 
estimation problems overestimation is more serious 
than underestimation, or vice-versa (Parsian & 
Kirmani (2002)). 

To deal with such cases, a useful and flexible class of 
asymmetric loss function (LINEX loss function (LLF)) 
was introduced by Varian (1975). The re- 
parameterized version of LLF (Singh et al. (2007)) for 
any parameter 0 is given as 

L(A) = e aA - a A - 1 ; 

a ^ 0 and A = (§ - ©)/© . (4.2) 

The sign and magnitude of ' a' represents the direction 
and degree of asymmetry respectively. The positive 
(negative) value of 'a' is used when overestimation is 
more (less) serious than underestimation. L(A) is 
approximately square error and almost symmetric if 
|a| near to zero. 

Bayes estimator of the parameter 0 under the LLF is 
denoted by 0 BL and obtained by simplifying the 
equation 
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Here the suffix p indicates the expectation taken 
under the posterior. Hence, the Bayes estimator under 
the LLF is 


0 B , — co 2 S . co 2 — 


1-exp 




n + d 


(4.3) 


JJ 


The expressions of the risk for the Bayes estimators are 
summarized in the following table under different loss 


Estimator 

Risk Under SELF 

0BS 

R (e Bs)=9 " + (<», -1 > 2 ) 

e BL 

R (0 bl)= e 2 (<n I + (co 2 -1 ) 2 ) 

Estimator 

Risk Under LLF 

0BS 

r(§bs)= acDj+a 1 

1-ato, 

e BL 

r(9bl)= a 0 2 + a 1 

l-aco 2 


The Minimax Bayes Risk Estimator 

The basic principle of Minimax Bayes risk estimation 
is to minimize the loss. Following Hodge and 
Lehmann (1950): 

Let CO = { F 0 : 0 e 0] be a family of distribution 
functions and C be a class of estimators of the 
parameter 0 . Suppose that c e C is a Bayes 
estimator against a prior distribution 71 (0) on the 

parameter space 0 . Then the Bayes estimator c is 
said to be the Minimax estimator if the risk function of 
the estimator C is independent on 0 . 

It is clear from the above table that the risks of the 
Bayes estimators under the SELF are not independent 
with the parameter 0 . Hence, the Bayes estimators 

/V /V 

0 BS and 0 BL are not the Minimax estimators. 

Therefore, under the SELF risk criterion the Minimax 
estimators do not exist. However, the risk of Bayes 

estimators 0 BS and 0 BL are independent of the 
parameter 0 under the LLF risk criterion. Hence, both 


the estimators 0 BS and 0 BL are Minimax estimators 
under the LLF loss criterion. 

The following statistical problem (Minimax Estimation) 
is equivalent to two person zero sum game between 
the Statistician (Player-II) and Nature (Player-I). Here 
the pure strategies of Nature are the different values of 
0 in the interval (0,co) and the mixed strategies of 
Nature are the prior densities of 0 in the interval 
(0,oo) . The pure strategies of Statistician are all 
possible decision functions in the interval (0,co). 

The expected value of the loss function is the risk 
function and it is the gain of the Player-I. Further, the 

Bayes risk of the Bayes estimator 0 B is defined as 

r(t,,0 b )=E 9 r( 0 b ) 

Here, the expectation has been taken under the prior 
density of the parameter 0 . If the loss function is 

continuous in both the estimator 0 B and the 

parameter 0, and convex in0 B for each value of 0 

* ^ * /v 

then their exist a measures T] and0 B for all 0 and0 B 
so that, the following relation holds: 

R ( n.9l)s R ( r ( n‘,e B )■ 

( * /v * \ 

T) ,0 B I is known as the value of the 

game, and T| and t) B are the corresponding 
optimum strategies of the Player I and II. In statistical 
terms T| is the least favourable prior density of 0 and 

the estimator 0 B is the Minimax estimator. In fact, the 
value of the game is the loss of the Player-II. Hence, 
the optimum strategy of Player-II and the value of 
game are given for present case as 


Optimum 

Strategy 

Corresponding 

Loss 

Value of 
Game 

^ BS 

LLF 

e~ a 

aco, + a - 1 

1-ato, 

9 BL = 

LLF 

e~ a 

aco., +a-l 

1-aa, 
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The Shrinkage Estimation 


Numerical Study 


Thompson (1968) suggested a procedure, which makes 
use of the prior information of the parameter 0 in 
form of a guessed value 0 O by shrinking the usual 
estimator towards the guess value of the parameter 
with the help of a shrinkage factor (see Prakash & 
Singh (2009)). 


The shrinkage procedure has been applied in 
numerous problems, including mean survival time in 
epidemiological studies (Harries & Shakarki, 1979), 
forecasting of the money supply (Tso, 1990), 
estimating mortality rates (Marshall, 1991), and 
improved estimation in sample surveys (Wooff, 1985). 
Following Thompson (1968), the shrinkage estimators 
based on the upper record values under the Bayesian 
Minimax criteria are defined as 

+(1 — (£> j ) 0 0 ;0 <COj<1. (6.1) 

and 

0 2 = co 2 0 BL +(l-co 2 )0 o ;0<co 2 <l. (6.2) 

It is observed later that for the considered set of 


parametric set of values, both COj and C0 7 lie between 
zero and one. Hence, they are chosen as the shrinkage 
factors. The risk of the shrinkage estimators 
0 j ; i = 1 , 2 under the LLF is given as 
o )_ e x P ((1 — COj) a 5 - a) 


R(e.l 


1 - a co , 2 


--aco ; 


-a8(l-co i ) + a-l ; i = 1, 2. 

Further, an unbiased estimator 0 LI for the parameter 

0 is also available. Therefore, a class of the unbiased 
estimator for the parameter 0 is constructed as 

T = k0 u ; keR + . (6.3) 

The value of k which minimizes the risk of the 
estimator T under the LLF is obtained as 

kitlin =— (l-e" a/2 ). 

a 

Hence, the improved class of estimator among the 
class (6.3) is 

t = k min 0 u , (6.4) 


with the risk under the LLF is given as 

R (t)= 2 (e~ a/2 -l)+a . 


Remark: 


The relative efficiency for the Bayes shrinkage 
Minimax estimator 0 l and 0 2 with respect to 
improved estimator T are defined as 



The expression of the relative efficiencies are the 
functions of a, d, 8 and 11 . For the selected set of 
values for a = 0.50, 1.00, 2.00; d = 0.50, 1, 2, 5, 10, 
50; 8 = 0.50(0.25) 1.50 and n = 05 (05) 20; the 

relative efficiencies have been calculated and 
presented here in the Tables 1-4. 

It is observed that both the Bayes Shrinkage Minimax 
estimators based on upper records value are uniformly 

performs better than the estimator T for all the 
considered set of parametric values. For both 
shrinkage estimators, the relative efficiency increases 
as sample size n increases for 8 < 1.00 and decreases 
otherwise. 

Similar trend has also been observed for ' d' . The 
relative efficiency attains maximum when 8 = 1 .00 
(except ford < 2.00) . Further, it is also noted that the 
relative efficiency decreases in the interval 8 > 0.75 
as the 'a' increases. 

Conclusions 

In present paper, we obtained the Shrinkage estimator 
based on the Bayes Minimax risk estimation criteria 
under the upper record values for the Rayleigh model. 
It is observed that on the basis of the relative efficiency, 
the proposed both Shrinkage Bayes Minimax 
estimators performs uniformly better than the 
improved estimator in a wide range of 8 which is 
defined here as the ratio between the true value and 
predicted (prior point) value of the unknown 
parameter 0 under the LLF criterion. Based on the 
magnitude of the relative efficiency, 0 , (Shrinkage 
Bayes Minimax estimator under LLF) is recommended 
for the interval 8 < 1 .00 and the estimator 0 , 
otherwise. 
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TABLE 3: RELATIVE EFFICIENCIES BETWEEN THE ESTIMATOR 0 ? AND T* FOR II =15 


d -» 

0.50 

1.00 

2.00 

5.00 

10.00 

25.00 

50.00 

'.50 

1.7784 

1.7844 

1.7956 

1.8233 

1.8560 

1.9077 

1.9423 

'.75 

5.3753 

5.4318 

5.5372 

5.8021 

6.1242 

6.6547 

7.0276 

.00 

98.464 

105.30 

119.66 

168.28 

267.74 

704.36 

1892.9 

.25 

16.878 

16.284 

15.282 

13.245 

11.422 

9.3649 

8.3499 

.50 

2.7707 

2.7194 

2.6301 

2.4355 

2.2440 

2.0030 

1.8721 

'.50 

1.7866 

1.7925 

1.8033 

1.8299 

1.8614 

1.9112 

1.9445 

'.75 

5.2161 

5.2696 

5.3694 

5.6201 

5.9248 

6.4264 

6.7787 

.00 

92.101 

98.458 

111.81 

156.98 

249.30 

654.12 

1754.9 

.25 

15.180 

14.638 

13.725 

11.870 

10.212 

8.3429 

7.4221 

.50 

2.3911 

2.3454 

2.2656 

2.0921 

1.9216 

1.7074 

1.5913 

'.50 

1.8011 

1.8065 

1.8166 

1.8415 

1.8709 

1.9173 

1.9483 

'.75 

4.9411 

4.9895 

5.0796 

5.3061 

5.5811 

6.0332 

6.3502 

.00 

81.340 

86.889 

98.537 

137.90 

218.21 

569.53 

1522.9 





